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Abstract. We have spectroscopically monitored the galactic Luminous Blue Variable HD 160529 and obtained an 
extensive high-resolution data set that covers the years 1991 to 2002. During this period, the star evolved from an 
extended photometric minimum phase towards a new visual maximum. In several observing seasons, we covered 
up to four months with almost daily spectra. Our spectra typically cover most of the visual spectral range with 
a high spectral resolution (A/ A A ~ 20 000 or more). This allows us to investigate the variability in many lines 
and on many time scales from days to years. We find a correlation between the photospheric Hel lines and the 
brightness of the star, both on a time scale of months and on a time scale of years. The short-term variations are 
smaller and do not follow the long-term trend, strongly suggesting different physical mechanisms. Metal lines also 
show both short-term and long-term variations in strength and also a long-term trend in radial velocity. Most 
of the line-profile variations can be attributed to changing strengths of lines. Propagating features in the line 
profiles are rarely observed. We find that the mass-loss rate of HD 160529 is almost independent of temperature, 
i.e. visual brightness. 
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1. Introduction 

A-type hypergiants and Luminous Blue Variables (LBVs), 
which often show spectra similar to A- hypergiants, are the 
visually brightest stars in galaxies. The early A-hypergiant 
HD 160529, which has long been known to exhibit the 



sorption lines of the order of 20 % and considerable pho- 
tometric variations were also found. The light variability 



was studied and described in detail by Sterken (1977) 
and [Sterken et al.| ( [l99lD . 



spectroscopic signatures of extreme luminosity (Merrill & 
Burwell 1943| ), has therefore long been considered as one 
of the visually brightest stars of the Galaxy, comparable 
to the brightest stars in external galaxies. 

The variability of HD 160529 was already established 
30 - 40 years ago through studies of its variable light and 
its spectrum at different epochs, which indicated changes 
in its s pectral type. The spectral variability of HD 160529 



HD 160529 was c lassified as a new galactic LBV by 
^terken et al. (1991) due to its brightness decrease of 
0.5 mag from 1983 to 1991 and an apparent spectral type 
change from A 9 to B 8. Apart from the long-term changes, 
the authors found pulsation-like variations in their photo- 
metric data with a quasi-period of 57 days and peak-to- 
peak amplitudes of 0.1 mag in b and y, while previous anal- 



was st udied extensively by Wolf et al. ( 1974), who found 
variations of the Balmer line profiles, radial velocity vari- 
ations with an amplitude of about 40 km s~^ and line 
splitting of some metal lines. Intensity variations of ab- 



yses suggested a possible 101.3 day period (Sterken 1981 
More recent photometry has been published by Manfroi^ 



eTaLl ( |l99ll , pMsj ) and pterken et al] (|l99|, p^95| ) 



From comparison with HD 269662 (=R110), a close 
photometric and spectroscopic counterpart of HD 160529 



O. 



Stahl, 



Send offprint requests to: 
. Stahl@lsw . uni-heidelberg . de 

* Based on observations collected at the European Southern 
Observatory at La Silla, Chile (proposals 69.D-0378, 269.D- 
5038) 



in the Large Magellanic Cloud, Sterken et al. (1991) de- 
rived an absolute visual magnitude of My = —8.9 and 
a distance of 2.5 kpc. They estimated the stellar parame- 
ters characterizing the phase of maximum visual bright- 
ness to be Teff X 8 000K, \ogg « 0.55, i?* 330 Rq 
and x 13 M©. The derived mass is in good agree- 
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ment with the mass determined from evolutionary tracks 
(Lamers et al. 1998). 

Sterken et ai.' ( |l99"ll) concluded that HD 160529 is lo- 
cated near the lower luminosity limit of the LBV insta- 
bility strip and possibly in an evolutionary phase after 
the Red Supergiant state, i.e. in a post-RSG evolutionary 
phase. 

HD 160529 is of much lower luminosity than more typ- 
ical LB Vs such as AG Car ( 3tahl et al. ^001 ) and shows 
much smaller variations of about 0.5 mag in the visual. 
According to van Gendcrcn ( 2001 ), it is at the low am- 
plitude limit of the "strong-active" SDor variables. The 
empirical relation between luminosity and amplitude of 
LBVs ( Wolf 1989| ) is in agreemen t with this. Accor d ing to 



the pei|lud-lumiiiusiLy iclaLiuii of Stnthuis Gten~(1995)' 



with its own cross-disperser, camera and CCD detector. 
The red channel is identical to the Flash instrument. The 
two spectral ranges of the Heros configuration cover the 
range from 3 450 A to 5 560 A in the blue channel and from 
5 820 A to 8 620 A in the red channel. The spectral reso- 
lution is about 20 000 for both Flash and Heros. The 
signal-to-noise ratio (S/N) of our spectra strongly depends 
on the wavelength and is lowest in the blue spectral range. 
For HD 160529 we typically used an exposure time of one 
to two hours for Flash and Heros with the ESO 50cm 
telescope. In good conditions, a S/N-ratio of at about 100 
is reached in the red spectral range. 

Further observations with significantly better S /N and 
spectral res olution were carried out with Feros ( Kaufei 



pcriud" uf wuU above 20 yi woulU Ue expecL 
is in rough agreement with the light curve published by 

pgil). 



juld be 



jLed. This 



et al.||2000D at the ESO 1.52m telescope between 1999 



Stcrkcn et al 



HD 160529 can be considered an intermediate case be- 
tween normal supergiants and LBVs. The closest coun- 
t erpart may be the LBV HD6884 (=R40) in the SMC 
( ^zeifert et a"Ll|l993|) or the LBV HD 269662 (=R110) in 
the LMC ( ^tahl et al.|p^ ). The variabihty of HD 160529 
is therefore of interest both in connection with typical 
LBVs such as AG Car ( ^tahl et al| |200l|) and with nor- 



mal late B - early A supergiants ( Kaufer et al. 19964|q) 



Detailed spectroscopic studies of LBVs, which cover 
both the short and and long time scales, are still very 



rare and only available for AG Car ( Stahl et al. 2001 ) . 
Studies of more objects with different physical parameters 
are needed in order to distinguish the general behaviour 
of LBVs from peculiarities of single objects. 

2. Observations 

2.1. Spectroscopy 

We have observed HD 160529 with several different echelle 
spectrographs and several telescopes from 1991, January 
to 2002, June. We have monitored HD 160529 with high 
time resolution over a time span of 2-4 months in four 
successive years, i.e. 1992-1995. In addition, a number of 
snapshot observations have been taken to extend the total 
time base covered. The epochs of observations and instru- 
ments used are summarized in Table ^. 

Most of the observations were obta ined wit h the fiber- 
linked echelle spectrographs Flash (|Mandel 



its modified version Heros (gtahl et al. 1996). These in 



19881 ) and 



struments were mainly used at the ESO 50cm telescope, 



but fo* 



few shorter runa also at the ESO 1.62m and 



2.2m t clcacopca. In the Flaoh configuration, a setting of 
the spectrograph was used that allowed the observation 
of the US and HeiA6678 lines with one exposure, i.e. the 
spectral range covered the wavelength range from about 
4 050 A to 6 7 80 A. The Flash data have been pubhshed 
on CD-ROM (|Stahl et al!||1995D. Heros is a modified ver- 



sion of Flash, where a beam splitter is used to divide 
the beam after the echelle grating into two channels, each 



and 2002. The Feros spectra obtained at the ESO 1.52m 
telescope cover the spectral range from about 3 600 A to 
9 200 A in 39 orders. The spectral resolution of Feros is 
about 48 000. The S/N depends strongly on wavelength 
and is highest in the red spectral region. The typical ex- 
posure times were about 10 minutes. 

In addition, a few spectra have been obtained with 
the Caspec echelle spectrograph at the 3.6m ESO tele- 
scope and the Ucles echelle spectrograph at the 3.9m 
AAT telescope. The Caspec spectra cover the spectral 
range from 3 800 A to 5 400 A with a spectral resolution of 
about 20 000. With Ucles, observations with three grat- 
ing positions were obtained and the resulting spectrum 
covers the total spectral range from 3 600 to 6 880 A. The 
spectral resolution is about 50 000. 

For all instruments, a built-in high-temperature incan- 
descent lamp and a Th-Ar lamp were used to obtain flat- 
field and wavelength-calibration exposures, respectively. 

In addition, we obtained on Sep. 24, 2001, low resolu- 
tion spectra (3A/pix in the wavelength range from 4 000 
to 8 000 A) with the Dfosc focal reducer spectrograph at 
the Danish 1.54m telescope at ESO, La Silla. In this case, a 
He-Ne lamp was used for the wavelength calibration. The 
Dfosc data are of much lower resolution than the echelle 
data and allow to study the equivalent width variations 
only. Since the spectra were obtained in non-photometric 
conditions, they can also not be flux-calibrated. 

The Ucles spectra were reduced with IRAF. All other 
spectra were reduced with ESO-Midas. The Caspec spec- 
tra were reduced with the standard echelle package of 
ESO-Midas ( |Ponz fc Brinksj |1983| ). For the Flash and 



Heros spectra a modified version of this package (Stahl 



et al. 1993) was used. The Feros package (Stahl et al 



1999), also running within ESO-Midas, was used for the 



reduction of the Feros spectra. The Dfosc spectra have 
been reduced with the long-slit package of ESO-Midas. 

In the years 1992-1995, the star was observed at least 
once per week and in some seasons up to once per night. 
In these runs, the sampling is sufficiently dense to cover 
variations with time scales of a few days. 
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Table 1. Summary of the spectroscopic observations of 
HD 160529. 



Instr. 


Telescope 


Sp. 


year 


Caspec 


ESO 3.6m 


2 


1991.08 


Flash 


ESO 50cm 


62 


1992.47-1992.75 


Flash 


ESO 50cm, 2.2m 


72 


1993.09-1994.41 


Flash 


ESO 50cm 


55 


1994.14-1994.49 


Heros 


ESO 50cm 


23b/27r 


1995.22-1995.42 


Ucles 


AAT 3.9m 


1 


1995.76 


Heros 


ESO 1.52m 


6 


1997.30-1995.32 


Feros 


ESO 1.52m 


2 


1999.53-1999.56 


Dfosc 


Danish 1.54m 


3 


2001.74 


Feros 


ESO 1.52rn 


2 


2002.49 52455 
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Fig. 1. Results of the photometry of HD 160529 in the 
Stromgren y band. The x-axis is in years. The epochs of 
spectroscopic observations are indicated by arrows, longer 
observing runs are indicated by double arrows. During our 
observations, the star evolved from photometric minimum 
to a new maximum. 

2.2. Photometry 

Differential photomet ry was carried out in the framework 
of the LTPV project (^terken| [l98^) using the ESO 50cm 
telescope (PMT photometry), the Dutch 90 cm telescope 
(CCD) and the Danish 1.54m telescope (CCD) at ESO, 
La Silla. Zero point corrections were applied but no proper 
standardization was attempted. The results based on the 
measurements in the Stromgren y band are displayed in 
Fig. I 



3. Variability 

3.1. Changes in spectral type 

The large photometric variations of LBVs are typically 
accompanied by changes of the spectral type, which can 



change from Of/WN- to A-hypergiant ( |Stah| |200l| ). For 
HD 160529, spectral types between B8Ia and A9Ia have 
been published in the literature. However, the spectrum of 
HD 160529 is peculiar and, since the spectral type depends 
on the exact classification criteria used, it is not clear if the 



differences in the published spectroscopic classifications 
are due to real variations. Most authors find a spectral 
type of A2~A3 ( |Wolf et al||l974^ , only two authors find 



a strongly deviating type: 



Houk (1982) finds A9Ia and 



^terken ct al.| (1991) infer spectral type variations from 
B8 to A9. 

It should be noted, t hat the spectral types published 
by 3terken et al. (1991) are not based on spectroscopic 



criteria, but they are inferred from two arguments: 

— The observed color variations variations indicate tem- 
perature changes from 8 000 - 10 000 K. This assumes 
that the theoretical colors apply to the extreme super- 
giant HD 160529. 

— The visual amplitude of 0.5 mag is due to a change 
in the bolometric correction at constant bolometric 
luminosity. This assumes that the bolometric correc- 
tions ( schmidt-Kaler 1982| ) for supergiants apply to 
HD 160529. 



The classification (A9Ia) of |Houk| ( |1982D is therefore 
the only spectroscopic classification s trongly devi ating 
from the A2-A3 range. The remark of Houk ( 1982 ): "H 
lines seem even weaker than in standards; partly filled in?" 
seems to indicate that this classification is partly based on 
the weakness of the Balmer lines, which are indeed influ- 
enced by emission. 

We therefore re-evaluated the evidence for strong 
changes in the spectral type of HD 160529. From our 
spectra we infer only minor changes in spectral type for 
the spectra obtained until 1999. For this period, we de- 
termined the spectral type of HD 160529 by comparing 
with the supergiant stars /3 0ri (B8Ia), HD 96919 (B9Ia), 
H D 92207 (AOIa ) and H D 100262 (A2Ia) from the sample 



of Kaufer ct al. (1997). The match of most lines is rea- 
sonably good with an approximate spectral type of A2Ia. 
The spectral type as determined from the strength of the 
metal lines is certainly always later than AO and earlier 
than A5Ia. Since most of the metal lines are forming partly 
in the expanding atmosphere of HD 160529, it is not clear 
how this spectral type relates to the star's effective tem- 
perature. This spectral type of A2Ia is in disagreement 
with the strength of HeiA5876, which would indicate a 
spectral type of B9Ia - A2Ia. The given range reflects the 
strong variations of this line (cf. Fig. H). 

We flnd that the spectrum before 1999 was signifi- 
cantly different from the maximum phase spectra of Woli 
et al. (1974). The spectral type was earlier (hotter) than 



in maximum phase. In addition, the emission lines are 
stronger in minimum phase. The temperature change is 
inferred mainly from the weakening of Tin and Cm lines 
in photometric minimum. However, since these lines form 
also in the stellar wind, the spectral type changes are 
unreliable indicators for the effective temperature. The 
HeiA5876 line may be a better temperature indicator than 
the metal lines, although the variable profile of HeiA5876 
indicates that this line is also not of purely photospheric 
origin. 
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Fig. 2. Radial velocity, FWHM and equivalent width of 
HeiA5876 as a function of time. The equivalent width of 
this line has been measured for other supergiants to 0.68 A 
for BSIa, 0.43A for B9Ia, 0.24A for AOIa and O.lOA for 
A2Ia. By this criterion, the star is coolest in 2002. 



Nevertheless, part of the spectral changes can certainly 
be attributed to changes in temperature. The strength of 
most metal lines strongly increased in 2002. In particular 
the CalA4227 line, which had equivalent widths around 
50 mA in earlier years, increased to about 270 mA in 2002 
(cf. Fig.| This line is a good temperature indicator for 
A stars. The strength of this line is slightly smaller than 



time was classified as about FOIa (Stahl ct al, 1990). The 
evidence is strong that the star in 2002 was considerably 
cooler than in the years before. The spectral type was 
earlier than FO, but much later than A2 and probably 
late A, confirming the spectral type variations reported 
earlier. 

3.2. Variations of line strengths and radial velocity 

Since HD 160529 developed a new S Dor phase during our 
observing campaign, strong spectroscopic changes are to 
be expected. 



HD 160529 




4224 4226 4228 4230 

Wavelength [A] 



Fig. 3. Variation of the CaiA4227 line. The strength of 
this line strongly increased in 2002. 



The general appearance of the spectrum of HD 160529 
changes very little within one observing season. Much 
stronger variations were observed during the whole cam- 
paign from 1992 to 2002. The lines of Cm and Tin are 
very temperature-sensitive in the temperature range of 
HD 160529 (i.e. about 8 000 to 12 000 K) and temperature 
changes should be reflected in the strength of these lines. 
Indeed, these lines exhibit particularly strong variations 
during our observing campaign. However, the line profiles 
are also highly variable, and is is therefore not clear 
that the changes in line strength are due to tem- 
perature changes only. Especially the short-term 
changes could be significantly affected by stellar 
wind effects (cf. e.g. Fig while the long-term 
changes, where these effects are averaged out, are 
more likely to be dominated by temperature vari- 
ations. Fig. 1^ shows an example of the spectral changes 
in Cm lines. The Cm are also sharper in 2002, which is 



also visible when comparing with Fig. 7 of Bterken ct al 



(1991): they show a spectrum obtained in 1990, where the 
Cm lines have a strong blue-shifted component. Compare 



with Fig. 4 of Wolf et al. (1974), where the Cm fines are 
always strong. Also by comp aring our equivale nt width 



measurements with Table 6 of |Wolf et al.| (1974), we find 
much lower values for Cm and Tin in most years. Only 



in the |LBV HD 269662 (=R110) in 1989, which at that ^Tal 



in 1999 and 2002, our values approach the values of |WoU 
( 1974 ) . The wea kening of Tin during minimum has 



also been reported by Sterken et al. (1991) 



Fig. ^ also shows strong variations of the H/3 
line, both in strength and line profile, between 
1992 and 2002. The line strength variations are 
probably caused by a density decrease in the stel- 
lar wind, caused by a radius increase and/or a de- 
crease of the stellar radius. The line profile shows 
that, while the terminal velocity is unchanged, the 
minimum of the absorption moved to smaller ex- 
pansion velocities. This could be due to decreased 
optical depth in the H/3 line or a decrease of the 
stellar wind velocity in the inner part of the wind. 



Otmar Stahl et al.: Spectroscopic monitoring of the LBV HD 160529 



5 




4860 
Wavelength [k] 



Fig. 4. The H/3 line and two Cm lines in 1992 (bottom) 
and 2002 (top). Note the much weaker emission and ab- 
sorption components of H/3 and the stronger metal ab- 
sorption lines in 2002 as compared with 1992. 



The HeiA5876 line is relatively symmetric and thus ap- 
pears well-suited to study the photospheric variations of 
HD 160529. The radial velocity of this line is highly vari- 
able, but also the FWHM and equivalent width of the line 
change. We fitted Gaussian profiles to all spectra and so 
obtained not only the radial velocity but also the residual 
intensity and line width as a function of time. The results 
are shown in Figs. |^ and |^. It is obvious from this figure 
that the line is not only variable in radial velocity but also 
in line width and intensity. Radial velocity and equivalent 
width are not correlated. The equivalent width and radial 
velocity are highly variable within each observing season, 
i.e. on time scales of a few months. In addition, the aver- 
age equivalent width seems to decrease from 1995 to 2002. 
This probably indicates a lower temperature in these years 
which is supported by e.g. Tin line measurements. The Ha 
strength also changed strongly in 1999-2002. The equiva- 
lent width was varying between about 14 and 17 A from 
1992 to 1997. In 1999 it was down to 10 A and further 
decreased to 7.5 A in 2001 and to 6.5 A in 2002. 

The relation of spectroscopic behaviour and light curve 
on short time scales can only be studied from the spectra 
obtained in 1992, when we have sufficient overlap in the 
spectroscopic and photometric data set (cf. Fig. ||). The 
sudden decrease of radial velocity seems to coincide with 
a peak of the visual brightness. The equivalent width of 
HeiA5876 and also the FWHM have a peak at the same 
time. Note that this correlation of line strength with vi- 
sual brightness on short time scales contrasts with the pro- 
nounced anti- correlation on long time scales (cf. Fig. 
These short-term variations of the HeiA5876 line 
could be due to pulsational variations or to stellar- 
wind effects. If they are due to pulsation, the vari- 
ations of the equivalent width imply an increase of 
the temperature with brightness. If we interpret 
the variations as changes of the stellar wind, the 
observations would suggest a strengthening of the 
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Fig. 5. Photometric and spectroscopic variations of the 
HeiA5876 line in 1992. Note the fast changes of almost 
15 km s~^ within only about 10 days and the relative 
constancy before and after this change. The equivalent- 
width curves seems to follow the light curve. 



wind with brightness, with a phase-delay of the 
radial velocity. 

Other lines which mainly originate in the photosphere 
are e.g. MgiiA4481 and the SiiiAA6347, 6371 lines. 

If variations in the strength of the Tin, Cm and Hei 
lines are due to photospheric temperature change of the 
star, we would expect the changes in Hei and the metal 
lines to be anti-correlated. Indeed, there is a weak evidence 
for this behaviour in some years, but the scatter is very 
large within individual years. In the long-term trend, the 
anti-correlation is clearly visible. 

This is probably due to the fact that the metal lines 
form mostly in the stellar wind (see Figs. ||and Only at 
the end of the observing run, the lines appear close to the 
systemic velocity and with a line width that is compatible 
with vsini. In other years, the radial velocity, line width 
and line profile variations indicate a stellar wind origin 
for these lines. Line splitting is quite common and has 



already been observed by Wolf et al. (1974). In particular 
in the first years of the observations, the radial velocities 
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Fig. 6. Radial velocity, FWHM and equivalent width of 
CrilA4824 as a function of time. The parameters have been 
determined by Gaussian fits to the line profiles, although 
the line profiles often are not Gaussian. Note the long- 
term trend in the radial velocity and the extremely strong 
changes in equivalent width on a time scale of weeks. The 
average equivalent width of about 0.3 A indicates a spec- 
tral type of about A2Ia, 0.1 A indicates AOIa. 



are more negative than observed by Wolf et al. ( 1974 ) in 
maximum visual brightness. 

3.3. Time series analysis 



Lamers et al. ( 1998 ) have previously analyzed the photo- 
metric data of HD 160529. For two seasons, they derived a 
period of 45 and 55 days, respectively. They explained the 
variations by non-radial g-mode pulsations of low order I 
and the period change by mode-switching. 

We examined the photometric data and radial velocity 
measurements for periodicities with the ESO-Midas TSA 



package using the method of Bcargle (1982). We limited 
the measurements to seasons with a significant number of 
measurements. For the photometry, these were the years 
1991-1993 and for the spectroscopy the years 1992-1995. 
We restricted the search to periods below 120 days, since 
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Fig. 7. Profile variations of CriiA4824 as a function of 
time in 1993. Note the strong change in line width and 
the occasional line-splitting. 



this is the approximate length of the observing runs. For 
the y-photometry we obtain the strongest peaks in the 
power spectrum at about 108, 83 and 68 days and for the 
radial velocity of Hei at 93 and 67 days. The equivalent 
width of this line shows two peaks at similar periods at 92 
and 67 days and in addition two peaks at 126 and 82 days. 
We checked that the results do not depend critically on the 
choice of the data set. If data sets from individual years 
are removed or added, the detected periods change, but we 
always find peaks at similar periods, i.e. vifithin a few 
days. This could indicate that the variations are 
not strictly periodic or not coherent over longer 
time scales. However, the small changes of the po- 
sition of the peaks of the power spectrum are more 
likely due to the irregular sampling of our data. 

It should also be mentioned that the radial velocity 
changes are definitely not sinusoidal. Although the time 
series analysis gives time scales of around 100 days, scales 
for the changes are of the order of 10 days. The time se- 
ries analysis measures the time scale on which such sudden 
changes of the radial velocity repeat. Examples of the ob- 
served variations are shown in Figs. |^ and ^. Similar events 
have been observed in other years as well. 

The similarity between the time scales of the photo- 
metric and spectroscopic variations suggests a physical 
connection. In 1992, when we obtained photometry and 
spectroscopy almost simultaneously, a peak in the light 
curve was observed to coincide with the a rapid change in 
radial velocity (Fig. |^). The equivalent width of HeiA5876 
also has a peak at the same time. In fact, the equivalent 
width seems to follow quite closely the visual light curve. 
Unfortunately, our observations do only allow to observe 
these relation in the 1992 observing seasons. Therefore, we 
can only conclude that the rapid changes in radial velocity 
may be related to the photometric variability. 
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3.4. Line profile variations 

Many lines in the spectrum of HD 160529 show peculiar 
line profile v ariations. While th is was already known from 
the work of Wolf et al. ( 1974 ), the time sampling of the 
spectra of these authors was not adequate to follow the 
changes in line profiles in detail. 

We searched in particular for propagating features in 
the line profiles, in order to study the stellar wind of 
HD 160529. Most of the lines are dominated by absorption 
and the observed variability seems mainly due to ad- 
ditional absorption. In order to increase the contrast, we 
therefore examined the residuals of the profiles with re- 
spect to a "maximum flux" profile, which was constructed 
from the time series by computing the highest profile for 
each wavelength bin in each observing season. This pro- 
cedure is very sensitive to noise. Therefore, a run- 
ning median was applied along the t-axis before 
computing the maximum flux. The subtraction of 
a "maximum flux" profile is more appropriate than 
subtracting a mean profile if the variations are due 
to additional absorption on top of a stable profile. 
The subtraction of a mean profile always intro- 
duces emission-like features in the residuals. 

A number of strong absorption lines do not show any 
emission component or a strong wind absorption 
and, therefore, superficially look like photospheric 
lines. An example of such metal lines arc the Sin lines 
at 6347, 6371 A. Profile variations as observed in 1993 are 
shown in Fig. ||. 

Most of the stronger metallic lines, in particular Fell 
lines, show pronounced P Cyg profiles. These lines show 
peculiar variations of the line profiles. An example of the 
variations is shown in Fig. Also in these profiles, most 
of the variability can be ascribed to additional absorption 
of variable strength, but with little change in radial veloc- 
ity. The variability is strongest at very low velocities 
around km s^. This component appears directly 
on top of the emission component. Therefore the 
P Cyg profile changes to a splitted absorption pro- 
file (see Fig. |l^, top) and later back to a P Cyg 
profile. Additional absorption at blue-shifted ve- 
locities of about 50 km s^ is also observed, but 
the features at low and high expansion velocities 
do not seem to be connected. This second, blue- 
shifted component seems to move to higher expan- 
sion velocities, but the changes are quite small. In 
some cases, the absorption features even seem to 
move towards longer wavelengths. Features which 
are clearly moving to higher expansion velocities, 
i.e. towards shorter wavelengths, are not often ob- 
served. Therefore, while we have some evidence for 
structures which originate in an accelerating stel- 
lar wind, the strongest variations occur close to the 
systemic velocity. A cool photospheric absorption 
layer may be responsible for these variations. 

The low-velocity feature is particularly strong in 1997, 
where it is clearly seen as a separate absorption feature. As 
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Fig. 8. Profile of the SiiiA6347 line in 1993. A "maximum 
flux" profile has been subtracted in the lower panel to 
enhance the contrast. 

an example we show in Fig. ^ the profile of the FeiiA4731 
line in the various observing seasons. In this figure it can 
be seen that the line was clearly blue-shifted from 1991 
to 1994 and later moved to longer wavelengths. The low- 
velocity feature could be related to the slightly red-shifted 



absorption features seen sometimes in LBVs (Wolf & Stah] 



1990; 3tahl 2001) 



A comparison of the line profile variations of lines 
which have a photospheric appearance (e.g. SiiiA6347 
in Fig. H) and lines which clearly show a P Cyg profile 
(e.g. FeiiA6248 in Fig. |l^) shows that the variations are 
very similar, although the mean profiles of SiiiA6347 and 
FeiiA6248 are very different. It should be mentioned, how- 
ever, that also the lines with pure absorption profiles show 
clear variations from year to year in radial velocity and line 
width. 

3.5. The resonance lines Nai, Can and Ki 

The NaiD lines show, apart from a strong interstellar com- 
ponent, clearly also a stellar contribution: a weak emission 
component and multiple variable absorption components. 
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Fig. 9. Profile of the FeiiA4731 line in various observing 
seasons from 1991 to 2002. The labels of the spectra denote 
the date in years. Note the strong low-velocity feature 

in 1997 and the strength and narrowness of the feature Fig. 10. Profile of the FeliA6248 line in 1993. A "maxi- 
in 1999. In 2002, the strength of the line further increases mum flux" profile has been subtracted in the lower panel 



and a strong blue wing develops. 



to enhance the contrast. Note the striking similarity with 
Fig. |. 



The variable absorption components are evidence for ac- 
celerating expanding shells. The maximum expansion ve- 
locity is typically around —150 km s^^ and reaches up 
to —200 km s^^. Discrete components are frequently ob- 
served in the velocity range from —60 to —150 km s^^. 

The changes in radial velocity of the components dur- 
ing one observing season are relatively small (up to 
about about 15 km s^, but typically below 10 km 
s~^) The changes from year to year are so large 
(of. Fig. that it is not clear hov^r the compo- 

nents are connected from year to year. Therefore, 
the development of the velocities on longer time- 
scales cannot reliably be determined from our ob- 
servations. Similar components may be present at 
lower velocities as well, but cannot be seen be- 
cause of the strong interstellar contribution, which 
saturates the absorption at velocities smaller than 
about —50 km s^^. Acceleration of the components in 
general is very small. In 1992, no significant acceleration 



was seen. In 1993 and 1995, an accelerating component 
was seen, most clearly in 1993 (cf. Fig. 12). Even in this 
year, we observed an acceleration of about 0.12 km 
s^^ d^^ at velocities from —75 to —90 km s^^. Another 
component at higher velocities (—115 km s~^) was stable 
in radial velocity. Both components first weakened and 
later strengthened again. In 1995, three components were 
observed (—125, —77 and —103 to —114) and one was 
accelerating. In 1994, two components were seen. Both 
components were very stable in radial velocity. The com- 
ponent at large expansion velocity (—120 km s^^) was 
weakening, while the component at smaller velocity (—75 
km s^^) was getting stronger with time. 

For CaiiH, K much less data is available. In particu- 
lar, these lines where not covered in the years 1992-1994, 
where we have obtained the most extended time series. 
The yearly mean spectra of the CaiiK line are shown in 
Fig. O. In cases where both the Can and the Nai lines 
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Fig. 11. Mean Nai profiles in the years 1992, 1993, 1994, 
1995, 1997, 1999 and 2002. The strongest variations are 
seen at the blue edge of the absorption component. The 
spectra are labeled with date in years. 
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Fig. 13. Mean Can profiles in the years 1991, 1995, 1997, 
1999 and 2002. The strongest variations are seen at the 
blue edge of the absorption component. The spectra are 
labeled with date in years. 
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pected acceleration is about 1 km s^^ d^^, i.e. significantly 
larger than what is observed. This implies that the line- 
splitting is probably not due to radially expanding shells 
moving with the stellar wind. This behaviour is similar 
to v^rhat is observed for discrete absorption compo- 



nents (DACs) in OB supergiants (of. e.g. FYillerton 
et al. 1997). The origin of these features and their 
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Fig. 12. Dynamic spectrum of NaiA5889 constructed from 
the spectra obtained in 1993. Note the accelerating feature 
around —75 to —90 km s^^. The acceleration of the feature 
is about 0.12 km s~^ d^^. 



are available, we find different velocities in both groups 
of lines. The blue-shifted components are much stronger 
in Can and therefore in most cases not as clearly sepa- 
rated. For the un-shifted component, which is mainly of 
interstellar origin, the Nai lines are stronger than Can. In 
Can, the interstellar contribution is split in two compo- 
nents at —26 and —4 km s~^. These components are only 
resolved in the Ucles and Feros spectra. This splitting 
is also marginally visible in the Nai lines. 

Overall, the picture is not clear. There is no dis- 
cernible long-term trend in the velocity or accel- 
eration of the components. From an expanding stellar 
wind {R = 200 R©, Woo = 200 km s'^, /3 = 4), the ex- 



slow acceleration is still not understood. 

The Ki line at 7699 A is another interstellar absorption 
line. This line has only been observed from 1995 to 2002. 
Due to the low abundance of K, the line is much weaker 
than Nai or Can. No blue-shifted components are seen. 
The line has two absorption components: a strong main 
component at —2 and a weaker component at —29 km 
s'^^. The components are clearly resolved only in the Feros 
spectra. From the two Feros spectra it appears that the 
fainter, blue component is slightly variable and therefore 
at least partly of stellar origin. 

4. Stellar parameters 

4.1. Distance 

The distance of HD 160529 is important for the deriva- 
tion of stellar parameters. Unfortunately, the distance is 
very uncertain: Sterken et al. (1991) estimate a distance 
of 2.5 kpc based on assigning to HD 160529 the same ab- 
solute visual magnitude as the LBV HD 269662 (=R110) 
in the LMC (My = —8.9). RllO is spectroscopically very 
similar to HD 160529 and also has a comparable photo- 
metric amplitude. In addition, an interstellar extinction 
of Ay — 3.4 mag (£'b-v = 1-1) is assumed. The empir- 
ical relation between luminosity and amplitude of LBVs 
( |Woq p^89| ) would give a lower bolometric luminosity of 
about A/boI = —8.3, if we assume an amplitude of 0.5 mag. 
This would also yield a slightly smaller distance estimate 
of about 1.9 kpc. If we assume that HD 160529 
has an absolute visual brightness of MboI — —9.5, 
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similar to HD 33579, the brightest A hypergiant 
of the Large Magellanic Cloud, the distance would 
be 3.5 kpc. 

A kinematic distance estimate is difficult for several 
reasons: the systemic velocity is difficult to determine be- 
cause of variable radial velocities and line profile varia- 
tions. In addition, the line-of-sight to HD 160529 is close 
to the galactic center and kinematic distances therefore 
have large errors. In the following, we will assume a dis- 
tance of 2.5 kpc but note that this distance has a rather 
large uncertainty of about 30%. 

4.2. Temperature changes 

The temperature of HD 160529 is another uncertain pa- 
rameter. Line strengths, derived e.g. from Sill, Hei and 
Fell or Till give contradictory results, when they are com- 
pared to normal supergiants or analyzed with static model 
atmospheres. This is certainly due to temperature strat- 
ification in the expanding atmosphere and stellar wind 
contributions in most absorption lines. 

The only metal which appears in different ionization 
stages is Mg. The Mgi/Mgii line ratio can be used to esti- 
mate the spectral type. The Mgi hues AA3832, 3838, 5173, 
5183 and the Mgii lines AA4481, 7877 were used. The re- 
sults are summarized in Table |^. For comparison, the AOIa 
star HD 92207, the A2Ia star HD 100262 and the F3Ia star 
HD 74180 were also measured. It should be noted that 
the line profiles of the measured lines also change with 
time. Therefore, the measured line ratios do not necessar- 
ily reflect changes in the effective temperature only. It is 
likely, however, that most of the variations of the 
line ratios reflect temperature changes. 

The near-IR also contains a large number of Ni lines, 
mainly in the wavelength range around 7 440 and 8 700 A. 
These lines are very strong, significantly stronger than 
in any of the late B/ear ly-A supergiants observed by 
Przybilla fc Butlei (2001). For the spectra taken between 
1995 and 1999, where the equivalent widths of these lines 
are relatively stable, we obtain for the lines AA7424, 7442, 
7468 mean equivalent widths of 325, 417 and 533 mA, re- 
spectively. In 2002 the lines are still stronger, with equiv- 
alent widths of 423, 529 and 614 mA, respectively. For 
HD 100262 (A2Ia), we measure 87, 141 and 196mA, re- 
spectively. The latter numbers are in very good agreement 



with the A2Ia star HD 111613 investigated by Przybilla & 
Butler] ( |2001| ). The much stronger hues for HD 160529 may 



indicate an enhanced abundance of N, but could also be 
due to non-LTE effects. 

Lines from neutral oxygen are strong in HD 160529. We 
measured an equivalent width of between 2.15 and 2.40 A 
for the near-IR triplet AA7771-5. For A stars, th e strength 
of this feat ure increases with decreasing \ogg ( Przybilla 
et al.|^OOOD . 

In principle, the temperature can also be derived from 
the continuum energy distribution. The optical energy dis- 
tribution is not very temperature sensitive, so the UV 



Table 2. Mgi and Mgii equivalent widths in mA. Note 
the strengthening of Mgi in 1997 and later. The line ra- 
tio Mgi/Mgii at most times gives a better match with 
HD 100262 (A2Ia) than with HD 92207 (AOIa). In 2002, 
the equivalent widths are exceptionally high and indicate 
a much later spectral type. The MgiA3832 line is blended 
in HD 74180. 



Mki 



Mgii 



season 


3832 


3838 


5173 


5183 


4481 


7877 


1992 






61 


61 


788 




1993 






64 


76 


783 




1994 






71 


66 


733 




1995 


156 


134 


89 


82 


761 


149 


1997 


164 


205 


137 


142 


706 


127 


1999 


164 


183 


116 


105 


714 


201 


2002 


397 


398 


378 


273 


834 


268 


HD 92207 


47 


43 


25 


29 


649 


104 


HD 100262 


161 


173 


124 


145 


667 


194 


HD 74180 




450 


553 


491 


766 


220 



energy distribution is needed. In the case of HD 160529 
this is difficult since the interstellar reddening is high and 
quite uncertain. However, low dispersion lUE spectra ob- 
tained in 1992 and 1979, i.e. during the visual minimum 
and maximum phase, respectively, show that the flux at 
about 2 500 A is lower by about a factor 2 in visual max- 
imum. This convincingly demonstrates that the tempera- 
ture was significantly lower at visual maximum. 

We tried to use synthetic spectra computed 
with the hydrostatic model atmosphere code 



Tlusty (Hubeny &: Lanz 1995; Lanz &: Hubeny 
2001| ) to determine the temperature of HD 160529 
and its changes. However, it turned out that the 
line spectrum of HD 160529 cannot be explained 
by hydrostatic models. Models with solar abun- 
dance which fit the HeiA5876 line in minimum re- 
quire about 11 000 K, while fits of metal lines such 
as, e.g., MgiiA4481 would require about 9 000 K. 
Increasing the He abundance reduces the discrep- 
ancy only slightly. We interpret this result as a 
strong indication of stratification effects in the ex- 
panding atmosphere. Non-LTE models which in- 
clude the atmospheric expansion and line blanket- 
ing are required. Wh ile such models are now avail- 
able for WR stars ( [Hillier Miller| |1998| ), they 
are so far not yet available for cooler stars like 
HD 160529. 

The changes of temperature and radius over the vari- 
ability cycle can be constrained by assuming that the bolo- 
metric magnitude of HD 160529 is constant. The change 
of 0.5 magnitudes in the visual brightness then reflects 
a change in the bolometric correction. Assuming a mini- 
mum temperature of 8 000 K and a corresponding radius 
of 330 R0, we derive a maximum temperature of 12 000 K 
and a corresponding radius of 150 R©. The change of the 
temperature is necessary for a corresponding change in 



Otmar Stahl et al.: Spectroscopic monitoring of the LBV HD 160529 



11 



bolometric correction. This conclusion follows assuming a 
black-body distribution or if wc use empirical bolomet- 
ric corrections (Schmidt-Kalei 1982). The large change in 
radius is required by the constancy of the bolometric lu- 
minosity. The temperatures correspond to spectral type 
of about B8 in visual minimum and A9 in visual maxi- 



mum ( schmidt-Kaler 198S), in agreement with the spec- 
troscopic evidence. 

4.3. Systemic velocity and rotational velocity 

The systemic velocity of HD 160529 is difficult to deter- 
mine since all absorption lines are influenced to some ex- 
tent by the stellar wind, their radial velocities are variable 
and different for different groups of lines. We therefore se- 
lected emission lines which appear little disturbed by ab- 
sorption to estimate the systemic velocity. The Fell lines of 
multiplets 40 and 46 appear well suited. From these lines 
we derive a radial velocity of —29 ± 3 km s~^. We also 
used the forbidden line NiiA6583 to estimate the systemic 
velocity. This line is flat-topped and probably forms at 
large distance from the star. The line center, as estimated 
from the bi-sector of the left and right edge, is at —8 km 
s~^, if we use a rest wavelength of 6583. 454A (Spyromilio 



1995). Unfortunately, the fine is faint and broad, so that 



the uncertainty of this measurement is large. 

The average velocity measured from HeiA5876 is about 
—23 km s~^ (see below), which is relatively close to this 
value. Other absorption lines such as the Cm multiplet 30, 
which appear as symmetric and narrow absorption lines 
in maximum, have velocities of around —20 km s~^. The 
red Ol and Nl lines give similar results: they vary between 
-27 and -19 km s'^ from 1995 to 2002. Because the 
radial velocity may be variable also for symmetric lines, 
it is not clear whether this velocity represents the true 
systemic velocity of HD 160529. 

Because of the variable line profiles, the rotational ve- 
locity vsini is also difficult to determine for HD 160529. 
Feros spectra taken in 1999 and 2002 show par- 
ticularly narrov^r absorption lines. We used these 
spectra to estimate an upper limit for usini. For 
this purpose, we fitted rotational profiles to some 
metal lines, taking into account the instrumental 
resolution of the spectrograph, but no other broad- 
ening mechanisms of the lines. In this way, we de- 
rive a vsini value of about 45 km s^^. This has to 
be considered as an upper limit, since stellar wind 
effects and other broadening mechanisms are most 
likely important. 

5. Modeling of Ha 

Mass-loss rates derived from radio observations are only 
little model-dependent. HD 160529 is one of the few stel- 
lar objects which are bright enough for current radio tele- 
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Fig. 14. Mean Ha profile in the years 1992, 1993, 1994, 
1995, 1997, 1999 and 2002. In the years 1992-1995, the 
strongest variations are seen at the blue edge of the ab- 
sorption component. After 1995, the emission decreased. 
The weakest emission is present in 2002, the second weak- 
est in 1999. 



wind at constant velocity. The authors adopt the spectro- 
scopic distance oi d = 2.5 kpc from Sterken et al. (1991) 
an electron temperature of Te = 5 000 K for the radio emit- 



ting region and a terminal velocity of Wqo = 180 km s to 
derive a mass-loss rate of Al = io~4-S'^^0-24 MQ-yr^'^. 

From our observations (cf. Fig. Q), we find that the 
emission-line flux of HD 160529 was slowly decreasing dur- 
ing our observing period, i.e. with increasing visual bright- 
ness. The variations during each season are much smaller 
than this long-term trend. 

In order to quantitatively investigate the mass-loss rate 
of HD 160529 versus time, we modeled the mean Ha line 
of the years 1992 and 2002, which show the strongest and 
weakest Ha emission, respectively. We used the multi- 



level line-transfer c ode described by Bastian (1982) and 
^tahl et al. (1983). A terminal velocity Voc of 200 km 
s~^ was used for both minimum and maximum. Note 
that stellar-wind theory predicts a decrease of Voc 
with decreasing surface gravity, i.e. with increas- 
ing visual brightness of LBVs. This is not easily 
observed in HD 160529, since there are no lines 
in the optical with a well-defined blue absorption 
edge. The absorption edge in Ha is variable, but 
this can also be explained by changes in the optical 
depth of the wind. A constant wind velocity can- 
not be excluded from our observations. We used a 
/3-type velocity law with (3—4. A velocity law with 



such a slow acceleration was found by Stahl et al 



scopes. Leitherer et al. (1995) observed the star at 8.64 and 
4.80 GHz and found that the measured flux is consistent 
with thermal emission from an optically thick expanding 



(2001) to fit the line profiles of AG Car in out- 
burst. The stellar radius and effective temperature were 
estimated from the photometry and the line spectrum of 
the given year. Then we only varied M and Vgys until a 
satisfactory fit was achieved. Solar abundances were as- 
sumed. Possible clumping was not taken into account. The 
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Table 3. Summary of the mass-loss rate determinations. 



= —10 km s and 



^sys 

both cases. 



= 200 km s where used in 



R/Rq 

Teff/K 

M/Moyr- 



150 330 
12 000 8 000 
710"" 110"^ 



results are summarized in Tab. js] and the fits are shown 
in Fig. 

Wc find that the mass-loss rate of HD 160529 is almost 
independent of temperature, i.e. visual brightness. Note 
that the derived mass-loss rate increases slightly from min- 
imum to maximum, although the equivalent width de- 
creases in this period. The reduced equivalent width is 
mainly due to the radius increase and ~ to a smaller ex- 
tent - to the decrease in temperature. The radius increase 
causes (at constant mass-loss rate) a density decrease 
which will decrease the equivalent width of the Balmer 
lines. Therefore, if the radius increase is smaller than de- 
rived from the assumption of constant bolometric lumi- 
nosity, the derived mass-loss rate could even slightly de- 
crease from minimum to maximum. The derived mass- 
loss rate is in good agreement with the results 
of Leitherer et al. ( [1995 ), which are based on ra- 
dio observations. We estimate that our mass-loss 
rates are accurate within about a factor of two. 
We therefore conclude that within the errors the 
mass-loss rate did not change significantly between 
minimum and maximum. 

A more advanced modeling of the stellar wind 
of HD 160529 w i th a code as e.g. described by 
Hillier Miller ( 1998 ) would certainly be worth- 
while. However, results of such models for cool 
stars such as HD 160529 are not yet published. 

6. Discussion and conclusions 

For the interpretation of the short-term variability of 
HD 160529, the relevant time scales are the rotational time 
scale Plot, the stellar- wind time scale Pwind and the pul- 
sation time scale Ppuis- 

The time scales are defined by the following relations: 

27ri?*/wbroak < -Prot < ^TrR^/v^ot (1) 

logPpuis = -0.275AfBoi - 3.918 log Teff + 14.543 (3) 

In Equation ^ the upper limit to Piot is estimated 
from vsini assuming sini = 1. This value can also be 
slightly larger because the derived v sin i value is an up- 
per limit. Note that it is expected that the rota- 
tional velocity is a function of phase. Because of 
the complex line profiles it could not be deter- 
mined in light minimum and is therefore set con- 
stant here. Because of the acceleration of the wind to its 
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Fig. 15. Fit of the Ha hne in 1992 and 2002. 

Table 4. Summary of time scale estimates (in days) for 
the minimum and maximum phase of HD 160529. 

min. (hot) max. (cool) 



Prot 59 ... 169 193 .. . 371 

Pwind >6 >13 

Ppuis 10 49 



terminal velocity in the extended stellar wind, the typical 
expansion time scale can be significantly longer than given 
by Equation ^. The pulsation period for the fundamental 
radial mode give n by Equat i on |^ is the empirical fitting 



formula given by iLovy et al| (|l984| ) 



We used the following parameters from Sterken et al 
( |1991D : Mboi = -8.9, = 13 Mq. In addition, wc used 
Voo = 200 km s^^ and Uj-ot = 45 km s^^ for both maximum 
and minimum phase. For the stellar radius and tempera- 
ture we adopted 150 Ro/12 000 K and 330 Rq/B 000 K for 
minimum and maximum phase, respectively. With these 
parameters, we obtain the time scale estimates given in 
Table |. 

Within the uncertainties, the time scales Pwind and 
Ppuis are compatible with the typical variation time scales 
of 50-100 days, while the rotational time scale is signif- 
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icantly longer. However, since rotational modulation 
can also lead to variations with an integer frac- 
tion of the rotational period, the time scale of the 
variations cannot be used to distinguish the pos- 
sible mechanisms of the short-term variations in 
HD 160529. The simultaneous variation of the optical 
brightness and the HeiA5876 line, however, is most easily 
explained by pulsations. Pulsations have also been sug- 
gested as a cause for the photometric micro-variations of 
LBVs (cf e.g. [Lamers ct aL||199^ ). 

The correlation of the equivalent width variations with 
visual brightness on short time scales is in marked contrast 
with the behaviour on long time scales. This strongly sug- 
gests that the physical mechanism for the variations on 
short time scales is different from the LBV-type varia- 
tions. 

It is rather surprising, that so far no strong dependence 
of the observed time scale on the visual brightness - and 
therefore the radius - has been detected. All plausible vari- 
ability mechanisms would predict s uch a dependence. It 
should be mentioned, however, that Bterken et al, (1991) 



derived a period of 57 days in minimum while the derived 
value in maximum was 101 days - about a factor of two 
longer. It is not clear, however, if the difference is due to 
a real period change. 

Most LBVs seem to increase their mass-loss rate in 
maximum. However, it is not clear that this is a gen- 



eral p roperty of LBVs (Leitherer 1997 ). Vink fc de Koter 
( 2002| ) have shown that the mass- loss rate of LBVs vcr 



sus temperature carries important diagnostic information. 
The mass-loss rate of AG Car versus phase has been stud- 
ied by [Stahl et al.| ( ^001] ). |Vink fc de Koteij (|00|) com- 
pared these results with their predictions for line-driven 
stellar winds and found a good agreement. They also pre- 
dict that the mass-loss rate of LBVs is a complicated func- 
tion of temperature, surface gravity and abundances. In 
particular, the mass-loss rate can even decrease with de- 
creasing temperature at temperatures below about 15 000 
K, in particular at low masses. This predicted behaviour 
nicely fits the observed behaviour of HD 160529, which is 
an LBV with very low temperature and mass. 

Although the mass- loss rates of LBVs can be predicted 
with some precision, the details of the mass-loss process 
are still not understood. In particular the complicated 
line-profile variations and the line splitting cannot be ex- 
plained by radially expanding features. Rotational modu- 
lation may be the cause for these variations. 
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